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ABSTRACT 



Observations show that the slope of the inner part of the Ha/ [Nil] rotation 
curves of disk galaxies is depressed by extinction: at fixed luminosity, the ob- 
served slope is in fact seen to depend on the disk inclination to the line of sight. 
Using a simple extinction model, we are able to reproduce well the observed 
trends. The model assumes an exponential distribution, both in the radial and z 
directions, identical for star-forming regions and dust. Fits to the data are op- 
timized by varying the scale height and scale length of absorbers, as well as the 
disk's central optical depth Tq, as seen face-on. The observations indicate that 
disk central opacity increases with total luminosity. Model fits for the most lumi- 
nous galaxies (brighter than M — 5 log/i = —21.4 in the I band) yield To ~ 3.5l2;o) 
near A = 0.66 /^m. The large uncertainty on the inferred value of To is due to the 
poorly known characteristics of the distribution of absorbers in the disk, as well 
as to the likelihood of strong departures from an exponential radial distribution 
near the galaxy centers. 

Subject headings: dust, extinction - galaxies; fundamental parameters - galaxies: 
halos -galaxies: spiral 
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1. 



Introduction 



Scattering and absorption by interstellar dust depresses the observed flux of galaxy 
disks at optical and infrared wavelengths. While the magnitude of the effect has been 
disputed, there is no doubt that the effect is strongly dependent on the inclination of a disk 
to the line of sight. A number of techniques have been used in observationally deriving 
extinction laws (see for example the volume edited by Davies and Burstein 1995, the review 
by Calzetti 2001 and Masters, Giovanelli & Haynes 2003 for details). Goad and Roberts 
(1981) first discussed a kinematical technique which we briefly describe as follows. 

Let V{r) be the rotational velocity of an axially symmetric disk at the distance r from 
its center and {x, y) be a set of Cartesian coordinates in the plane of the disk. If the disk 
is thin and its symmetry axis is inclined by an angle i to the line of sight, the observed 
component of velocity which intercepts the disk at (x, y) is 



where V^c accounts for non-circular motions and the x axis is oriented along the disk's 
apparent major axis. The observed rotation curve, as derived for example from a long-slit 
Ha spectrum positioned along the major axis, will depart from the prescription given by 
Eqn. 1 as due to seeing, instrumental resolution, averaging across the slit width, the finite 
thickness of the disk and extinction occurring within the disk itself. As a realistically 
thick disk approaches the edge-on perspective, lines-of-sight along the major axis sample 
regions of increasingly broad range in y, yielding a velocity distribution with a peak 
velocity contributed by parcels of gas asymmetrically distributed with respect to y = 0: if 
extinction is important, only foreground parts of the disk contribute to the emission and 
the factor (x/ a/x^ + y'^) < 1 depresses the velocity distribution observed aX r = x. Goad 
& Roberts noted how this tapering effect may, in opaque edge-on disks, yield observed 
rotation curves resembling solid-body behavior, independently of the true shape of V{r). 



V\\ = V{r) 



X 



sin i + 



(1) 
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Bosma et al. (1992) tested this technique on two edge-on systems, NGC 100 and NGC 891. 
They compared HI synthesis and Ha observations and concluded that the disk of NGC 
100 is transparent, while in the case of NGC 891 they could not exclude the possibility of 
extinction in the inner parts of the disk. Their conclusions were however strongly influenced 
by the limited spatial resolution of the HI data. Prada et al. (1994) compared long-slit 
spectra of NGC 2146 in the optical and near IR and reported evidence for extinction in 
the inner parts of the galaxy. More recently, Giovanelli & Haynes (2002, hereafter GH02) 
applied the same technique in a statistical fashion, detecting clear evidence for the effect 
of extinction at 0.66 /^m. Using a sample of more than 2000 Ha rotation curves, they 
found that the opacity in the inner disks of spiral galaxies is luminosity dependent, in a 
manner previously found by purely photometric means (Giovanelli et al. 1995; TuUy et al. 
1998). The observational data of GH02 are corroborated by the detailed radiative transfer 
calculations of Baes et al.(2004), which also obtain dramatic increases in the apparent inner 
scale length of rotations curves as disks approach the edge-on perspective. 

In this report, we use a simple model for the dust gas distribution in spiral disks, in 
order to reproduce the findings of GH02 and in so doing to obtain quantitative inferences 
on the characteristics of the disk opacity. In Section 2 we summarize the observational 
results that are to be modelled. In Section 3 we describe the disk model to be fitted to the 
data, while the best fit parameters are discussed in Section 4. In Section 5 we present our 
conclusions. Throughout this work, distance dependent quantities are scaled according to a 
Hubble parameter Ho — lOO/i km s~^ Mpc~^ . 

2. Rotation Curve Shape and Disk Inclination 

GH02 fitted Ha/ [Nil] rotation curves of spiral galaxies with the parametric model 

Fpe = K(l - e-'-/'^-)(l + ar/r^e) (2) 
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where K scales the amphtude of the rotation curve, Vpe yields a scale length for the inner 
steep rise, and a corresponds to the slope at large r. In their Table 1, for different luminosity 
classes as measured in the I-band they list mean values of a, hvp^ and of the ratio between 
hvpe and the scale length of the disk hght for relatively face-on systems. In their Figure 1, 
for each luminosity class the dependence of hvpe on disk inclination i is shown. For a given 
luminosity class, with increasing i, < hvpe > remains at first constant; as i approaches 70°, 
< hrpe > starts increasing, the rate of increase being larger for more luminous systems. In 
the latter, < hrpe > high inclination has values more than double the values observed at 
low inclination. The inclination dependence was interpreted as due to the increase in the 
opacity of spiral disks with inclination, while the luminosity dependence was attributed to a 
combination of increased disk size and higher interstellar medium metallicities in the more 
luminous systems. By contrast, the outer slope < o; > is unaffected by disk inclination. 



3. Disk Extinction Model 

As done by other authors who have studied disk extinction models (e.g (Byun et al. 
1994; Xilouris et al. 1999; Misiriotis et al. 2000), we assume that the emissivity per unit 
volume is produced in the disk by a distribution of isotropic emitters which is exponential 
both in the radial direction r in the plane of the disk, and in the direction z perpendicular 

to it, so that 

j(r, z) = jo e-'/''' e-\'^l'^ (3) 

where Vd and Zd are respectively the scale length and the scale height, and we ignore the 
effects of dumpiness and azimuthal variations. We will need to distinguish between scale 
lengths and scale heights of stars, dust and Ha-emitting regions, respectively referred to 
by subscripts "*" , "d" and "a" . We will assume that the dust distribution in the disk is 
exponential both in r and z^ with scale length and scale height of the dust as those of the 
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HcK light at A = 0.66 //m. That assumption greatly simplifies the analysis; it is justified since 
the light distribution in our case is that of the Ha produced in star forming regions, which 
is likely to mimick that of dust clouds. It is however a simplistic assumption, which cannot 
fully reproduce the more complex characteristics of the distribution of star-forming regions 
in galaxies. We use it here because of its simplicity, and discuss later some of its possible 
limitations. The absorption coefficient per unit volume for an exponential distribution can 
thus be written as 

K{r, z) = Ko e-''/'"''-'' e-l^l/^'*''^ (4) 

The optical depth at the disk center, when seen face-on, is then Tq = 2«;o^d,d, while 
the optical depth as a function of the spatial coordinate along the line of sight s is 
t{s) = K,{s')ds'. For a randomly oriented disk, the optical depth to any parcel of 
emitting gas at location s is fully specified once To, ra^d, Zd,d and i are given. The disk's 
intrinsic axial ratio for each component, go = Zd/vd, is related to the apparent ratio between 
minor and major radii b/a and the inclination to the line of sight via 

cos'i^[(b/af-ql]/(l-q'o), (5) 

which is needed to relate the model to the observations, for GH02 show the inclination 
dependence of the rotation curve shape by plotting hvpe vs. logiQ^a/b). GH02 derived axial 
ratios from I-band surface photometry, a band for which the light originates in late disk 
stars; the go,* of that light is very likely to be significantly larger than that of the dust and 
of the star-forming regions. In Giovanelli et al. (1994,1995) they derived a value of go,* of 
0.13 for the I-band. The largest data sets of Ha surface photometry arc those of Ryder & 
Dopita (1994) and Koopmann, Kenney & Young (2001). Ryder and Dopita have compared 
stellar and Ha scale lengths, showing the latter to be significantly larger than the former; 
from their figures 4 and 5, we infer rd,a = (1.9 ± 0.4:)rd,*, where the stellar scale length 
corresponds to the I-band. 
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The contribution to the observed flux along a given line of sight by a parcel of gas 
at location s is proportional to j{s)e~'^^^\ The Ha line profile along a given line of sight 



velocity, given a rotation curve V{r). In analogy with the measurement practice by which 
the rotation curve of a disk is extracted from a slit positioned along the disk major axis, 
in our model we assign to each line of sight the velocity corresponding to the peak of the 
simulated Ha line profile. 

We use the model as follows: 

• Given a disk galaxy luminosity class, we obtain an "unextincted" rotation curve V{r) 
of the type described by Eqn. 2, using the average values of the structural parameters 
as given in Table 1 of GH02. 

• We select a set of disk parameters ra^d, qo,d — the same as for the Ha light — and 
To and build an "observed" rotation curve, that takes into account the effects of 
extinction as produced by the adopted disk model. We then measure hrpg in the thus 
modified rotation curve. 

• By progressively changing the inclination angle of the disk, we obtain a curve that 
describes the change of hrpe with cos(i), which we compare with the observed data 
as shown in Figure 1 of GH02. In order to show the dependence with cos(i) in the 
observed data, we use the Eqn. 5. 

• The comparison between model and data is optimized via the minimization of 



where [rpe,i, cos(i)i] is the ith point of the data to be fitted, and Uj is its associated 
uncertainty. The procedure yields a best fit of unconstrained values for ra^d, g^d and 



through the disk will be oc j{s)e '^^^^ds, where s is mapped onto the line of sight 
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To, or best fit values of any of those parameters, when the others are constrained. 
Since GH02 express the disk inchnation in terms of the axial ratio of the starlight in 
the I-band, we use Eqn. (5) with = Qo,* = 0.13 to obtain cos(i). 

4. Model Fitting Results 

The GH02 analysis of rotation curves separates the available sample in six luminosity 
classes. The statistically more significant trend of h rpg with inclination takes place for the 
two brightest luminosity classes: Class I includes galaxies with M — 5 log h < —22 and Class 
II includes galaxies with —22.0 < M — blogh < —21 A. The trends in hvpe with inclination 
are shown by the data points in Figure 3. The average stellar scale lengths in the I band are 
respectively hvd,* — 4.3 kpc and hrd,* — 3.2 kpc, for Classes I and II, as given by Haynes et 
al. 1999. Here, we restrict our analysis to those two classes. For the less luminous classes, 
extinction is small, the inclination dependence of hvpe is marginal or nonexistent and the 
value of To that would be derived is unreliable. 

Figures 1 and 2 show simulated rotation curves, as obtained using our disk extinction 
model, Eqn. 2 and the parameters hsted in Table 1 of GH02, for the two luminosity classes 
of interest. In Figure 1, the "unextincted" rotation curve for Lumnosity Class I, normalized 
by a factor VoSini is shown as a thick line in each of the two panels. In panel (a), the 
effect of changing the disk inclination is shown, with thin lines exhibiting rotation curves 
for values of i between 70° and 86°, Ai = 4°, at a fixed value of Tq = 3.0, for hrd^ — 8-6 
kpc and go,d = 0.05, plotted versus the radial distance in units of the stellar scale length 
Td,*- In panel (b), for the same values of hrd,d and g^d, the effect of varying the central 
disk optical depth is shown, whereby for a fixed inclination oi i — 86°, the values for 
To = 0.5, 1.0, 2.0, 5.0, 10.0, are shown in the corresponding curve. Figure 2 displays the 
analogous model curves for Luminosity Class II, for which we have assumed hrd,d = 6.4 kpc 
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and go,d = 0.05. In all panels, the intersection of the horizontal dashed line and the rotation 
curve identifies the corresponding value of Tpe for that curve. Note how even for disks of 
moderate opacity, the value of the rotational velocity within a few scale lengths from the 
center can be drastically altered by extinction, in the case of highly inclined and luminous 
systems. 

Figure 3 shows the GH02 data on the variation of hrpe with inclination, and the 
corresponding model trends as derived from our disk extinction model, for varying values of 
To (1 to 10, Ato = 1). In each case, g^d and hrd,d are the same values used in Figures 1 and 
2. The models plotted with a thick line correspond to Tq — 5.0. 

In order to verify the variance in the best fit value Tq produced by different values of 
the assumed go,d, Figure 4 shows the data for the two luminosity classes and, separately, 
simulated hrpg vs. cos(i) relations for Tq between 1 and 10 and for g^d = 0.03 (panels a, d), 
0.05 (panels b, e) and 0.07 (panels c, f). Similarly, Figure 5 shows the variance introduced 
by different assumptions on the dust disk scale length hrd,d- for qo,d = 0.05, the model 
curves were computed for hrd,d — 7.0 kpc (panel a), 8.0 kpc (panel b), and 9.0 kpc (panel 
c) for Luminosity Class I, and 5.5 kpc (panel d), 6.5 kpc (panel e) and 7.5 kpc (panel 
f) for Luminosity Class II. For given Tq, disk inclination and r^ ^i, a higher value of go,d 
implies that extinction along a given path length involves a wider interval of disk radii, thus 
resulting in a more "smeared" rotation curve (i.e. longer hrpe). Similarly, for given To, disk 
inclination and go,d, a longer dust scale length also implies that extinction along a given 
path length involves a wider interval of disk radii, thus a longer hvpg. Both effects are seen 
in the simulations. 

Unconstrained fits of the type described by Equation (6) can yield simultaneously 
values for To, qo,d and rd,d- However, values of qo,d and rd,d are not independent: lengthening 
the disk will obviously reduce its axial ratio. We can thus maintain < rd,d > and/or < go,d > 



-10- 



fixed and obtained constrained fit values for < To > for each Luminosity Class. What is a 
reasonable range of values for < q^^^ > and < ra^d > ? 

While photometric measurements provide a fair assessment of the mean and variance 
in the value of hrd,* for a given luminosity class (the scatter about the observed mean 
values is approximately 0.08 in logr^,*), the distribution in disk intrinsic axial ratios is quite 
uncertain, for g^* is a difficult quantity to measure and even its mean value for a given type 
of galaxy relies on poorly constrained statistical inferences from the observed distribution of 
axial ratios and from a few observations of edge-on systems. The mean value of go,* = 0.13 
was adopted by Giovanelli et al. (1994) for the I band. The extreme disk population of star 
forming regions and dust should be significantly thinner than that of the solar-type disk 
stars contributing most of the light at I band. Ryder & Dopita (1994) give a comparison 
between ra,* and rd,a-i which we assumed to be the same as rd,d- No equally accessible 
comparison between go,* and go,ci is available. In Table 1, we give estimates of To obtained 
assuming the tabulated values for r^^d and go,^, based on the Ryder & Dopita (1994) data 
and heuristic assumptions on the dust layer intrinsic axial ratio. The error assigned to this 
determination is based on a combination of (a) the standard error propagation analysis 
derived from the impact of the GH02 data quality, (b) uncertainty in the conversion of Vd,* 
to Td^d, estimated to be 20% on the basis of the Ryder & Dopita data, and (c) uncertainty 
in the value of go,d, assumed to hover between 0.035 and 0.06 for both Luminosity Classes. 
The largest source of uncertainty is clearly that deriving from go^^. An additional source of 
uncertainty is introduced by the often-seen departures from an exponential form, in the 
radial distribution of Ha especially near the centers of galaxies: a cursory inspection of the 
Ryder & Dopita and the Koopmann et al. samples indicates that about half of galaxies 
exhibit a cusp in the Ha distribution, within the inner 1-2 kpc. This region coincides 
with that within which r^e is measured. Thus, values of To inferred with the assumption 
that the radial distribution of Ha, and thus of dust, is exponential throughout the disk 
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may underestimate the opacity of the inner regions of the disk and overestimate it in the 
outer regions. While it is clear that extinction plays a role in affecting the shape of the 
inner parts of Ho; rotation curves of bright spirals, the technique outlined in this paper 
for its measurement requires higher quality estimates of the distribution of absorbers than 
currently available, for To to be pinned down to better than within a factor of 2 or so. 

5. Conclusions 

A simple model, whereby dust and HIT regions are distributed exponentially in both 
the radial and z directions, with equal scale height and scale length, can reproduce the 
observed variation in the inner slope of Ha/ [Nil] rotation curves of disk galaxies. The 
effect is due to internal extinction. For reasonable values of the ratio between scale height 
and scale length, the model can constrain central disk opacities for the most luminous 
galaxies, in which the effect of extinction is the strongest. However, pinning down the 
exact value of central opticall depth of the disk at 0.66 //m is difficult, due to uncertainty 
in the form of the distribution of absorbers. If Ha and the absorbers have scale length 
about 1.9 times that of the I band starlight and scale height to scale length ratio of ~ 0.05, 
we infer central optical depth of disk, in the face-on perspective of To — 3.5, for galaxies 
brighter than M — 51og/i = —21.4 in the I band. However, the uncertainty in the radial 



Table 1. Inferred Disk Central Optical Depth To 



Lum Class 


/iTd,* (kpc) hrd,d (kpc) 




To 


I 

II 


4.3 8.2 ±1.7 
3.2 6.4 ±1.2 


0.045 ± 0.015 
0.050 ±0.015 


9 c:+2.5 

q p;+4.0 
0.0_2.o 
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and ^-distributions of absorbers constrains the inferred values of To to an uncertainty no 
better than a factor of 2. For less luminous spiral galaxies, extinction becomes progressively 
less important, and the model cannot usefully constrain the disk opacity 

Support for this work was provided by the National Science Foundation through grant 
number GF- 1003-00 from the Association of Universities for Research in Astronomy, Inc., 
under NSF Coorporative Agreement No. AST-9613615. 
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Fig. 1. — Simulated rotation curves for luminosity class I. The solid line corresponds to the 
enextincted disk. Panel (a) has rotation curves for fixed Tq = 3 and varying disk inclination, 
from i — 70° to i — 86°. In panel (b), rotation curves are for fixed i — 86° and varying To, 
as indicated. The intersection of the horizontal, dashed line and each rotation curve occurs 
near the radius hr„e. 
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Fig. 2. — Simulated rotation curves for luminosity class II. The solid line corresponds to the 
enextincted disk. Panel (a) has rotation curves for fixed Tq = 3 and varying disk inclination, 
from i — 70° to i — 86°. In panel (b), rotation curves are for fixed i — 86° and varying To, as 
indicated in figure 1. The intersection of the horizontal, dashed line and each rotation curve 
occurs near the radius /ir„p. 
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Fig. 3. — GH02 data on the variation of hvpe with inchnation for luminosity classes I and 
II (data points) and simulated inclination dependences as obtained by our model, varying 
values of Tq between 1 and 10 (bottom to top solid lines). The thick line corresponds to 
To = 5 in each panel. 
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Fig. 4. — hrpe vs. cos{i) for Luminosity Class I (Luminosity Class II) varying To between 
1 and 10 (bottom to top solid lines), for panel a(d) go,d = 0.03, panel b(e) go,d = 0.05 and 
panel c(f) qo,d = 0.07. The thick line corresponds to To = 5 in each panel. 
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Fig. 5. — hrpe vs. cos(z) for Luminosity Class I varying To between 1 and 10 (bottom to 
top solid lines), for (a) hr^^d = 7.0 kpc, (b) hr^^ci = 8.0 kpc and (c) /ir^^^d = 9.0 kpc, and 
Luminosity Class II for (d) hr^^d = 5.5 kpc, (e) /ir^d = 6.5 kpc and (f) /ir^d = 7.5 kpc. The 
thick line corresponds to Tq = 5 in each panel. 



